Big Bang nucleosynthesis with a non-Maxwellian distribution 
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The abundances of light elements based on the big bang nucleosynthesis model are calculated using 
the Tsallis non-extensive statistics. The impact of the variation of the non-extensive parameter q 
from the unity value is compared to observations and to the abundance yields from the standard big 
bang model. We find large differences between the reaction rates and the abundance of light elements 
calculated with the extensive and the non-extensive statistics. We found that the observations are 
consistent with a non-extensive parameter q = llo.ilj indicating that a large deviation from the 
Boltzmann-Gibbs statistics (q = 1) is highly unlikely. 



I. INTRODUCTION 

The cosmological big bang model is in agreement with 
many observations relevant for our understanding of the 
universe. However, comparison of calculations based on 
the model with observations is not straightforward be- 
cause the data are subject to poorly known evolutionary 
effects and systematic errors. Nonetheless, the model is 
believed to be the only probe of physics in the early uni- 
verse during the interval from 3 — 20 min, after which the 
temperature and density of the universe fell below that 
which is required for nuclear fusion and prevented ele- 
ments heavier than beryllium from forming. The model 
is inline with the cosmic microwave background (CMB) 
radiation temperature of 2.275 K p], and provides guid- 
ance to other areas of science, such as nuclear and parti- 
cle physics. Big bang model calculations are also consis- 
tent with the number of light neutrino families N v = 3. 
According to the numerous literature on the subject, 
the big bang model can accommodate values between 
N, y = 1.8 — 3.9 (see, e.g., Ref. From the measure- 

ment of the Zq width by LEP experiments at CERN one 
knows that N v = 2.9840 ± 0.0082 [3]. 

In the big bang model nearly all neutrons end up in 
4 He, so that the relative abundance of 4 He depends on 
the number of neutrino families and also on the neutron 
lifetime r„. The sensitivity to the neutron lifetime af- 
fects Big Bang Nucleosynthesis (BBN) in two ways. The 
neutron lifetime r„ influences the weak reaction rates be- 
cause of the relation between r„ and the weak coupling 
constant. A shorter (longer) r„ means that the reaction 
rates remain greater (smaller) than the Hubble expansion 
rate until a lower (larger) freeze-out temperature, having 
a strong impact on the equilibrium neutron-to-proton ra- 
tio at freeze-out. This n/p-ratio is approximately given 
in thermal equilibrium by n/p = exp[— Am/ksT] ~ 1/6, 
where fcs is the Boltzmann constant, T the temperature 
at weak freeze-out, and Am is the neutron-proton mass 
difference. The other influence of r„ is due to their de- 
cay in the interval between weak freeze-out (t ~ 1 s) 
and when nucleosynthesis starts (t ~ 200 s), reducing 
the n/p ratio to n/p ~ 1/7. A shorter r n implies lower 



the predicted BBN helium abundance. In this work we 
will use the value of r n = 878.5 ± 0.7 ± 0.3 s, according 
to the most recent experiments [4] (a recent review on 
the neutron lifetime is found in Ref. |5J. Recently, the 
implications of a change in the neutron lifetime on BBN 
predictions have been assessed in Ref. [5]. 

The baryonic density of the universe deduced from the 
observations of the anisotropies of the CMB radiation, 
constrains the value of the number of baryons per photon, 
n, which remains constant during the expansion of the 
universe. Big bang model calculations are compatible 
with the experimentally deduced value of from WMAP 
observations, n = 6.16 ± 0.15 x 10~ 10 [7]. 

Of our interest in this work is the abundances of light 
elements in big bang nucleosynthesis. At the very early 
stages (first 20 min) of the universe evolution, when 
it was dense and hot enough for nuclear reactions to 
take place, the temperature of the primordial plasma de- 
creased from a few MeV down to about 10 keV, light 
nuclides as 2 H, 3 He, 4 He and, to a smaller extent, 7 Li 
were produced via a network of nuclear processes, re- 
sulting into abundances for these species which can be 
determined with several observational techniques and in 
different astrophysical environments. Apparent discrep- 
ancies for the Li abundances in metal poor stars, as mea- 
sured observationally and as inferred by WMAP, have 
promoted a wealth of new inquiries on BBN and on stellar 
mixing processes destroying Li, whose results are not yet 
final. Further studies of light-element abundances in low 
metallicity stars and extragalactic H II regions, as well as 
better estimates from BBN models are required to tackle 
this issue, integrating high resolution spectroscopic stud- 
ies of stellar and interstellar matter with nucleosynthesis 
models and nuclear physics experiments and theories [8] . 

The Maxwell-Boltzmann distribution of the kinetic en- 
ergy of the ions in a plasma is one of the basic inputs for 
the calculation of nuclear reaction rates during the BBN. 
The distribution is based on several assumptions inherent 
to the Boltzmann-Gibbs statistics: (a) the collision time 
is much smaller than the mean time between collisions, 
(b) the interaction is local, (c) the velocities of two par- 
ticles at the same point are not correlated, and (d) that 



2 



energy is locally conserved when using only the degrees of 
freedom of the colliding particles (no significant amount 
of energy is transferred to and from collective variables 
and fields). If (a) and (b) are not valid the resulting effec- 
tive two-body interaction is non-local and depends on the 
momentum and energy of the particles. Even when the 
one-particle energy distribution is Maxwellian, additional 
assumptions about correlations between particles are nec- 
essary to deduce that the relative- velocity distribution is 
also Maxwellian. Although the Boltzmann-Gibbs (BG) 
description of statistical mechanics is well established in 
a seemingly infinite number of situations, in recent years 
an increasing theoretical effort has concentrated on the 
development of alternative approaches to statistical me- 
chanics which includes the BG statistics as a special limit 
of a more general theory [5j (see also, [TO])- Such theories 
aim to describe systems with long range interactions and 
with memory effects (or non-ergodic systems). A very 
popular alternative to the BG statistics was proposed 
by C. Tsallis OUT], herewith denoted as non-extensive 
statistics (for more details on this subject, see the exten- 
sive reviews [TTHT5] ). Statistical mechanics assumes that 
energy is an "extensive" variable, meaning that the total 
energy of the system is proportional to the system size; 
similarly the entropy is also supposed to be extensive. 
This might be justified due to the short-range nature 
of the interactions which hold matter together. But if 
one deals with long-range interactions, most prominently 
gravity; one can then find that entropy is not extensive 

In classical statistics, to calculate the average values 
of some quantities, such as the energy of the system, the 
number of molecules, the volume it occupies, etc, one 
searches for the probability distribution which maximizes 
the entropy, subject to the constraint that it gives the 
right average values of those quantities. As mentioned 
above, Tsallis proposed to replace the usual (BG) en- 
tropy with a new, non-extensive quantity, now commonly 
called the Tsallis entropy, and maximize that, subject to 
the same usual constraints. There is actually a whole infi- 
nite family of Tsallis entropies, indexed by a real-valued 
parameter q, which quantifies the degree of departure 
from extensivity (one gets the usual entropy back again 
when q = 1). It was shown in many circumstances that 
the classical results of statistical mechanics can be trans- 
lated into the new theory [13] . The importance of these 
families of entropies is that, when applied to ordinary 
statistical mechanics, they give rise to probabilities that 
follow power laws instead of the exponential laws of the 
standard case (for details on this see [13 ). In most cases 
that Tsallis formalism is adopted, e.g. Ref. [JJ5], the 
non-extensive parameter q is taken to be constant and 
close to the value for which ordinary statistical mechan- 
ics is obtained (q = 1). Some works have also probed 
sizable deviations of the non-extensive parameter q from 
the unity to explain a variety of phenomena in several 
areas of science [11] . 

In the next sections, we shown that the Maxwell- 



Boltzmann distribution, a cornerstone of the big bang 
and stellar evolution nucleosynthesis, is strongly mod- 
ified by the non-extensive statistics if q strongly devi- 
ates from the unity. As a consequence, it also affects 
strongly the predictions of the BBN. There is no "a pri- 
ori" justification for a large deviation of q from the unity 
value during the BBN epoch. In particular, as radiation 
is assumed to be in equilibrium with matter during the 
BBN, a modification of the Maxwell distribution of ve- 
locities would also impact the Planck distribution of pho- 
tons. Recent studies on the temperature fluctuations of 
the cosmic microwave background (CMB) radiation have 
shown that a modified Planck distribution based on Tsal- 
lis statistics adequately describes the CMB temperature 
fluctuations measured by WMAP with q = 1.045±0.005, 
which is close to unity but not quite |20| . Perhaps more 
importantly, Gaussian temperature distributions based 
on the BG statistics, corresponding to the q — »• 1 limit, 
do not properly represent the CMB temperature fluctu- 
ations [ID]. Such fluctuations, allowing for even larger 
variations of q might occur during the BBN epoch, also 
leading to a change of the exponentially decaying tail of 
the Maxwell velocity distribution. 

Based on the successes of the big bang model, it is fair 
to assume that it can set strong constraints on the limits 
of the parameter q used in a non-extensive statistics de- 
scription of the Maxwell-Boltzmann velocity distribution. 
In the literature, attempts to solve the lithium problem 
has assumed all sorts of "new physics" [8]. The present 
work adds to the list of new attempts, although our re- 
sults imply a much wider impact on BBN as expected 
for the solution of the lithium problem. If the Tsallis 
statistics appropriately describes the deviations of tails 
of statistical distributions, then the BNN would effec- 
tively probe such tails. The Gamow window (see figure 
[T]) contains a small fraction of the total area under the 
velocity distribution. Thus, only a few particles in the 
tail of the distributions contribute to the fusion rates. 
In fact, the possibility of a deviation of the Maxwellian 
distribution and implications of the modification of the 
Maxwellian distribution tail for nuclear burning in stars 
have already been explored in the past [2"TM24j . As we 
show in the next sections, a strong deviation from q = 1 
is very unlikely for the BNN predictions, based on com- 
parison with observations. Moreover, if q deviates from 
the unity value, the lithium problem gets even worse. 



II. MAXWELLIAN AND NON-MAXWELLIN 
DISTRIBUTIONS 

Nuclear reaction rates in the BBN and in stelar evolu- 
tion are strongly dependent of the particle velocity dis- 
tributions. The fusion reaction rates for nuclear species 
1 and 2 is given by (av) 12, i.e., an average of the fusion 
cross section of 1+2 with their relative velocity, described 
by a velocity distribution. It is thus worthwhile to study 
the modifications of the stellar reaction rates due to the 
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modifications introduced by the non-extensive statistics. 

A. Non-extensive Statistics 

Statistical systems in equilibrium are described by the 
Boltzmann-Gibbs entropy, 



Sbg = -kB^pilnpi, 



(1) 



where ks is the Boltzmann constant, and pi is the proba- 
bility of the i-th microstate. For two independent systems 
A, B, the probability of the system A + B being in a state 
i + j, with i a microstate of A and j a microstate of B, 
is 



A+B A B 
Pi+j Pi P 3 ■ 



(2) 



Therefore, the Boltzmann-Gibbs entropy satisfies the re- 
lation 



Sa+b — Sa + Se 



(3) 



Thus, the entropy based on the Boltzmann-Gibbs statis- 
tic is an extensive quantity. 

In the non-extensive statistics [5] , one replaces the tra- 
ditional entropy by the following one: 
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(4) 



where q is a real number. For q = 1, S q = Sbg- Thus, 
the Tsallis statistics is a natural generalization of the 
Boltzmann-Gibbs entropy. 
Now it follows that 



S q (A + B) = S q (A)+S q (B) 



(1-9) 



S q (A)S q (B). (5) 



The variable q is a measure of the non-extensivity. Tsal- 
lis has shown that a formalism of statistical mechanics 
can be consistently developed in terms of this general- 
ized entropy [13] . 

A consequence of the non- extensive formalism is that 
the distribution function which maximizes S q is non- 



Maxwellian [25H27] . For q = 1, the Maxwell distribu- 
tion function is reproduced. But for q < 1, high energy 
states are more probable than in the extensive case. On 
the other hand, for q > 1 high energy states are less 
probable than in the extensive case, and there is a cutoff 
beyond which no states exist. 



where a is the fusion cross section, v is the relative veloc- 
ity of the y-pair, iVj is the number of nuclei of species i, 
fi is the reduced mass of i+j, T is the temperature, S(E) 
is the astrophysical S-factor, and r\ = ZiZje 2 /hv is the 
Sommerfeld parameter, with Zi the i-th nuclide charge 
and E = /j,v 2 /2 is the relative energy of i + j. 

The energy dependence of the reaction cross sections 
is usually expressed in terms of the equation 



a(E)= S ^cxp[-27rr 1 (E)}. 



We write 2-kt\ = b/y/E, where 



b = 0.9898ZiZ,-/A MeV 1/2 , 



(7) 



(8) 



where A is the reduced mass in amu. The factor 1 + 8ij 
in the denominator of Eq. ^ corrects for the double- 
counting when i = j. The S-factor has a relatively weak 
dependence on the energy E, except when it is close to a 
resonance, where it is strongly peaked. 



C. Non-Maxwellian Distribution 

The non-extensive description of the Maxwell- 
Boltzmann distribution corresponds to the substitution 



f(E) -> f q (E), where [13] 



f q {E) 



9->l 



exp 



q -^E 
k B T 

E 



< E < oo. 



(9) 



If q - 1 < 0, Eq. (|9| is real for any value of E > 0. 
However, if q — 1 > 0, f(E) is real only if the quantity in 
square brackets is positive. This means that 



< E < 



k B T 

9-1' 
< E, 



if q > 1 
if q < 1. 



(10) 



Thus, in the interval < q < 1 one has < E < oo and 
for 1 < q < oo one has < E < E max = k B T/(q - 1). 



B. Maxwellian Distribution 

In stars, the thermonuclear reaction rate with a 
Maxwellian distribution is given by [25] 



1 + 5, 



1 + 5, 



dES(E) exp 



13 
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knT 



k B T 



(6) 



With this new statistics, the reaction rate becomes 

N,N, 



,J ~ ! • -i, , '''- 



and the rate integral, I q , is given by 

pE moI 



dES(E)M q (E, T), 



(11) 



(12) 
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Figure 1: Modified Gamow distributions M q (E,T) of 
deuterons relevant for the reaction 2 H(d,p) 3 H at Tg = 1. The 
solid line, for q — 1, corresponds to the use of a Maxwell- 
Boltzmann distribution. Also shown are the results when us- 
ing non-extensive distributions for q — 0.5 (dotted line) and 
q = 2 (dashed line). 



where the "modified" Gamow energy distribution is 



M q (E,T) = A(q,T)[l 



-E 



-b/VE 



x exp 



-0.98982^ 



0.08617T 9 

A 

E 



(13) 



is the non-extensive Maxwell velocity distribution, 
E m ax — oo for < q < 1 and E max = k B T/(l - q) 
for 1 < q < oo, and E in MeV units. A(q,T) is a nor- 
malization constant which depends on the temperature 
and the non-extensive parameter q. 



D. Non-Maxwellian distribution for relative 
velocities 

It is worthwhile to notice that if the one-particle energy 
distribution is Maxwellian, it does not necessarily imply 
that the relative velocity distribution is also Maxwellian. 
Additional assumptions about correlations between par- 
ticles are necessary to deduce that the relative-velocity 
distribution, which is the relevant quantity for rate cal- 
culations, is also Maxwellian. This has been discussed in 

where non-Maxwellian distri- 



details in Refs. EEJ EI 
butions, such as in Eq. (13), were shown to arise from 
non-extensive statistics. 

Here we show that, if the non-Maxwellian particle ve- 
locity distribution is given by Eq. ([9]), then a two-particle 
relative can be modified to account for the center of mass 
recoil. Calling the kinetic energy of a particle £7,-, this 



distribution is given by 
f(i) - ( i _ q ~ 1 E- 
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The two-particle energy distribution is fq f q 
now exponentiate the Tsallis distribution. 
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We 

(15) 



and the product fq 12 ^ = fq 1 ^ reduces to 



f(12) 
J q 



exp 



1,1 ! ' - ^ 



1 - ^E 2 
kT ' 



(16) 

Since Ei — rmvf/2, and thus, E\ + E 2 — ^v 2 /2 + 
MV 2 /2, where ^ is the reduced mass of the two particles, 
M = mi +TO2, v is the relative velocity, and V the center 
of mass velocity, the product inside the natural logarithm 
can be reduced to 
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(17) 



Thus, the two-body distribution factorizes into a prod- 
uct of relative and center of mass parts 



4 12 ) (v, V, T) = ( v , T)& m ~> (V, T) 



(18) 



where 



fjf el \v,T) = A rel (q,T) \l- 



1 — q \xv 
kT ~Y 

1 - q MV 2 \ T 
~kf 2~ 



(19) 



with the normalization constants obtained from the con- 
dition, 



d 3 vd 3 Vf^. 12) (v,V,T) = 1. 



(20) 



Because the distribution factorizes, the unit normaliza- 
tion can be achieved for the relative and cm. distribu- 
tions separately. The distribution needed in the reaction 
rate formula is, therefore, 

f q (v,T) = [ d 3 Vf( 12 Hv,V,T) = & el) (v,T), (21) 



which attains the same form for as the absolute distribu- 
tion. 

In the limit q — > 1 the two-particle distribution reduces 
to a Gaussian, with the last term in the left-hand-side of 



Eq. (17) dropping out, 



jf 2 ) = .4(g, T) exp 
as expected. 



liv 2 /2 + MV 2 /2 
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Figure 2: The relative difference ratio (n. 
non-extensive, n 9 , and extensive, n 
curves are for Fermi-Dirac statistics, n + , and dashed curves 
are for Bose-Einstein statistics, n~. For both distributions, 
we use /i = 0. Results are shown for q — 2 and g = 0.5, with 
T g = 10. 



E. Equilibrium with electrons, photons and 
neutrinos 

One of the important questions regarding a plasma 
with particles (i.e., nuclei) described by the non-extensive 
statistics is how to generalize Fermi-Dirac, Bosc-Einstcin, 
and Tsallis statistics, to become more unified statistics 
with the distribution for the particles. This has been 
studied in Ref. [30] . where it was shown that a simi- 
lar non-extensive statistic for the distribution can be ob- 
tained for fermions and bosons is given by 



nf(E) 



1 



1 - (q - 1 



kT 



(23) 



± 1 



where /i is the chemical potential. This reproduces 
the Fermi distribution, n + , for q — >• 1 and the Bose- 
Einstein distribution for photons, n - , for [i = and 
q 1. Planck's law for the distribution of radiation is 
obtained by multiplying nT in Eq. ( 23 ) (with fj, = 0) 
by TaJ 2 1 (Att 2 c 2 ) , where E = %u). The number den- 



sity of electrons can be obtained from n + in Eq. (23) 
with the proper phase-factors depending if the electrons 
are non-relativistic or relativistic. Normalization factors 
A ± {q,T) also need to be introduced, as before. 



the nuclear reactions during the BBN. In fact, the en- 
hancement of the nuclear reaction rates due to electron 
screening were shown to be very small [3T] . The electron 
Fermi energy for these densities are also much smaller 
than kT for the energy relevant for BBN, so that one can 
also use /i = in Eq. (23) for n + . 

In figure[2]we plot the the relative difference ratio (n J — 
n )/n between non-extensive, n q , and extensive, n = 
n 9 _>.i, statistics. For both distributions, we use /i = 0. 
The solid curves are for Fermi-Dirac (FD) statistics, n + , 
and dashed curves are for Bose-Einstein (BE) statistics, 
n~ . We show results for q = 2 and q = 0.5, with Tg = 
10. One sees that the non-extensive distributions are 
enhanced for q > 1 and suppressed for q < 1 , as compared 
to the respective FD and BE quantum distributions. The 
deviations for the FD and BE statistics also grow larger 
with the energy. For example, the non-extensive electron 
distribution is roughly a factor 2 larger than the usual 
FD distribution at E e = 1 MeV, at T 9 = 10. 

While we don't obtain numerical results with modified 
Fermi-Dirac and Bose-Einstein distributions here, it can 
be expected that these generalizations will have a strong 
influence on the freezout temperature and the neutron 
to proton, n/p, ratio. A numerical study of this prob- 
lem may be presented in another paper. The freezout 
temperature, occurs when the rate, T ~ (<Jv), for weak 
reaction v e + n — > p + e~ becomes slower than the ex- 
pansion rate of the Universe. Because during the BBN 
the densities of all particles, including neutrinos, are low 
compared to kT, the chemical potential /i can be set to 
zero in the calculation of all reaction rates. The adop- 
tion of non-extensive quantum distributions such as in 
Eq. ( 23 ) will lead to the same powers of the tempera- 



ture as those predicted by the FD distribution and the 
Bose-Einstein distribution. For example, Planck law for 
the total blackbody is U oc T 4 , being form invariant with 
respect to non-extensivity entropic index q which deter- 
mines the the degree of non-extensivity |32j . This result 
means that the weak decay reaction rates do not depend 
on the non-extensive parameter q. The freezout temper- 
ature and n/p ratio remain the same as before. 

More detailed studies have indicated that Planck's law 
of blackbody radiation and other thermodynamical quan- 
tities arising from non-extensive quantum statistics can 
yield different powers of temperature than for the non- 
extensive case [33-36 . If that is the case, then a study 
of the influence of non-extensive statistics on the weak- 
decay rates and electromagnetic processes during BBN is 
worth pursuing. 



The electron density during the early universe varies 
strongly with the temperature. At T 9 = 10 the electron 
density is about 10 32 /cm 3 , much larger than the elec- 
tron number density at the center of the sun, nl un ~ 
10 26 /cm 3 . The large electron density is due to the e + e~ 
production by the abundant photons during the BBN. 
However, the large electron densities do not influence 



F. Thermodynamical equilibrium 

The physical appeal for non-extensivity is the role 
of long-range interactions, which also implies non- 
equilibrium. Accepting non-extensive entropy means 
abandoning the most important concept of thermody- 
namics, namely the tendency of any system to reach 
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equilibrium. This also means that the concept of non- 
extensivity means renouncing to the second law of ther- 
modynamics altogether! 

The comments above, which seem to be shared by part 
of the community (see, e.g., p£fl are worrisome 

when one has to consider a medium composed of parti- 
cles obeying classical and quantum statistics. It is not 
clear for example if the non-extensive parameter q has 
to be the same for all particle distributions, both classi- 
cal and quantum. Even worse is the possibility that the 
temperatures are not the same for the different particle 
systems in the plasma. 

In the present work, we will avoid a longer discussion 
on the validity of the Tsallis statistics for a plasma such 
as that existing during the BBN. We will only consider 
the effect of its use for calculating nuclear reaction rates 
in the plasma, assuming that it can be described by a 
classical distribution of velocities. This study will allow 
us to constrain the non-extensive parameter q based on 
a comparison with observations. 
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Figure 3: S-factor for the reaction 2 H(d,p) 3 H as a function of 
the relative energy E and of the temperature Tg. The data 
are from Refs. [42H46] . The solid curve is a polynomial fit to 
the experimental data. 



III. REACTION RATES DURING BIG BANG 
NUCLEOSYNTHESIS 

Based on the abundant literature on non-extensive 
statistics (see, e.g., |2"TH2~4] ). we do not expect that 
the non-extensive parameter q differs appreciably from 
the unity value. However, in order to study the influence 
of a non-Maxwellian distribution on BBN we will explore 
values of q rather different than the unity, namely, q = 0.5 
and q = 2. This will allow us to pursue a better un- 
derstanding of the nature of the physics involved in the 
departure from the BG statistics. In figure [I] we plot the 
Gamow energy distributions of deuterons relevant for the 
reaction 2 H(d,p) 3 H at T 9 = 1. The solid line, for q = 1, 
corresponds to the use of the Maxwell-Boltzmann distri- 
bution. Also shown are results for non-extensive distri- 
butions for q = 0.5 (dotted line) and q — 2 (dashed line). 
One observes that for q < 1, higher kinetic energies are 
more accessible than in the extensive case (q = 1). For 
q > 1 high energies are less probable than in the exten- 
sive case, and there is a cutoff beyond which no kinetic 
energy is reached. In the example shown in the figure for 
q = 2, the cutoff occurs at 0.086 MeV, or 86 keV. 



We will explore the modifications of the BBN elemen- 
tal abundances due to a variation of the non-extensive 
statistics parameter q. We will express our reaction rates 
in the form Na{o"v) (in units of cm 3 mol" 1 s" 1 ), where 
Na is the Avogadro number and (av) involves the in- 
tegral in Eq. ^ with the Maxwell distribution f(E) 
replaced by Eq. First we show how the reaction 

rates are modified for q =/= 1. 

In figure [3] we show the S-factor for the reaction 
2 H(d,p) 3 H as a function of the relative energy E. Also 
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Figure 4: Reaction rates for 2 H(d,p) 3 H as a function of the 
temperature Tg for different values of the non-extensive pa- 
rameter q. The rates are given in terms of the natural loga- 
rithm of Na(cv) (in units of cm 3 mol -1 s _1 ). Results with 
the use of non-extensive distributions for q — 0.5 (dotted line) 
and q = 2 (dashed line) are shown. 

shown is the dependence on Tg (temperature in units of 
10 9 K) for the effective Gamow energy 

E = E = 0.l22{ZfZ]A) 1/3 Tg /3 MeV, (24) 

where A is the reduced mass in amu. The data are from 
Refs. 42 - JS] and the solid curve is a chi-square polyno- 
mial function fit to the data. 



Using the chi-square polynomial fit obtained to fit the 
data presented in figure [3j we show in figure [4] the reac- 
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tion rates for 2 H(d,p) 3 H as a function of the tempera- 
ture Tg for two different values of the non-extensive pa- 



rameter q. The integrals in equation (12 1 are performed 



numerically. For charge particles, a good accuracy (witi- 
hin 0.1%) is reached using the integration limits between 
E - 5AE and E + 5AE, where AE is given by Eq. ( [25] ) 
below. The rates are expressed in terms of the natural 
logarithm of Na(ctv) (in units of cm 3 mol -1 s _1 ). The 
solid curve corresponds to the usual Maxwell-Boltzmann 
distribution, i.e., q = 1. The dashed and dotted curves 
are obtained for q = 2 and q = 0.5, respectively. In both 
cases, we see deviations from the Maxwellian rate. For 
q > 1 the deviations are rather large and the tendency 
is an overall suppression of the reaction rates, specially 
at low temperatures. This effect arises from the non- 
Maxwellian energy cutoff which for this reaction occurs 
at O.O86T9 MeV and which prevents a great number of 
reactions to occur at higher energies. 



> 

CD 



CO 




E (MeV) 



Figure 5: S-factor for the reaction 7 Li(p,a) 4 He as a function 
of the relative energy E and of Tg. The data are from Refs. 
[47H57] . The solid curve is a chi-square function fit to the 
data using a sum of polynomials plus Breit-Wigner functions. 



For q < 1 the nearly similar result as with the Maxwell- 
Boltzmann distribution is due to a competition between 
suppression in reaction rates at low energies and their 
enhancement at high energies. The relevant range of en- 
ergies is set by the Gamow energy which for a Maxwellian 
distribution is given by Eq. ( 24 ) and by the energy win- 
dow, 



AE = 0.2368(Z? Z]A) 1/g Tq /6 MeV, 



(25) 



5/6 



which for the reaction 2 H(d,p) 3 H amounts to 0.2368T 
MeV. This explains why, at Tg = 1, the range of rele- 
vant energies for the calculation of the reaction rate is 



shown by the solid curve in figure [1] For q < 1 the 
Gamow window AE is larger and there is as much a con- 
tribution from the suppression of reaction rates at low 
energies compared to the Maxwell-Boltzmann distribu- 
tion, as there is a corresponding enhancement at higher 
energies. This explains the nearly equal results shown in 
figure [4] for q = 1 and q < 1. This finding applies to 
all charged particle reaction rates, except for those when 
the S-factor has a strong dependence on energy at, and 
around, E = Eq. But no such behavior exists for the 
most important charged induced reactions in the BBN 
(neutron- induced reactions will be discussed separately) . 
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Figure 6: Reaction rates for 7 Li(p,a) 4 He as a function of the 
temperature Tg for two different values of the non-extensive 
parameter q. The rates are given in terms of the natural 
logarithm of Na{o~v) (in units of cm 3 mol" 1 s" 1 ). Results 
with the use of non-extensive distributions for q — 0.5 (dotted 
line) and q = 2 (dashed line) are shown. 



The findings described above for the reaction 
2 H(d,p) 3 H arc not specific but apply to all charged par- 
ticles of relevance to the BBN. We demonstrate this with 
one more example: the 7 Li(p,a) 4 He reaction, responsi- 
ble for 7 Li destruction. In figure [5] we show the S-factor 
for this reaction as a function of the relative energy E. 
One sees prominent resonances at higher energies. Also 
shown in the figure is the dependence of the reaction on 
T 9 . The data are from Refs. [4TH5T] and the solid curve 
is a chi-square function fit to the data using a sum of 
polynomials plus Breit-Wigner functions. 

Using the chi-square function fit obtained to fit the 
data presented in figure [5j we show in figure [6] the re- 
action rates for 7 Li(p,a)^Ie as a function of the tem- 
perature Tg for two different values of the non-extensive 
parameter q. The rates are given in terms of the natural 
logarithm of Na(<tv) (in units of cm 3 mol" 1 s" 1 ). The 
solid curve corresponds to the usual Maxwell-Boltzmann 
distribution, i.e., q — 1. The dashed and dotted curves 
are obtained for q — 2 and q = 0.5, respectively. As with 
the reaction presented in figure [4j in both cases we see 
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deviations from the Maxwellian rate. But, as before, for 
q = 2 the deviations are larger and the tendency is a 
strong suppression of the reaction rates as the temper- 
ature decreases. It is interesting to note that the non- 
Maxwellian rates for q = 0.5 are more sensitive to the 
resonances than for q > 1. This is because, as seen in 
figure [T] for q < 1 the velocity distribution is spread 
to considerably larger values of energies, being therefore 
more sensitive to the location of high energy resonances. 
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Figure 7: Spectral function M q (E, T) for protons and neu- 
trons relevant for the reaction p(n,7)d at Tg — 0.1 (upper 
panel) and Tg = 10 (lower panel). The solid line, for q = 1, 
corresponds to the usual Boltzmann distribution. Also shown 
are non-extensive distributions for q — 0.5 (dotted line) and 
q = 2 (dashed line). 



We now turn to neutron induced reactions, which are 
only a few cases of high relevance for the BBN, notably 
the p(n,7)d, 3 He(n,p)t, and 7 Be(n,p) 7 Li reactions. For 
neutron induced reactions, the cross section at low ener- 
gies is usually proportional to where v — \JlmE jfi 
is the neutron velocity. Thus, it is sometimes appropriate 
to rewrite Eq. ^ as 

S(E) _ R(E) 



a(E) = 



E 



4e 



(26) 



where R{E) is a slowly varying function of energy simi- 
lar to an S-factor. The distribution function within the 



reaction rate integral (12 1 is also rewritten as 



M q (E,T) = A(q,T)f q (E) = A(q,T)(l- 



g-1 



E 



(27) 



The absence of the tunneling factor exp(— b/v E) in Eq. 
( 27 ) inhibts the dependence of the reaction rates on the 
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Figure 8: The energy dependence of R(E) = S(E)\ /r E for the 
reaction 7 Be(n,p) 7 Li. is shown in figure Isl The experimental 
data were collected from Refs. [58H62] . The solid curve is a 
function fit to the experimental data using a set of polynomi- 
als and Breit-Wigner functions. 



In figure [7] we plot the kinetic energy distributions of 
nucleons relevant for the reaction p(n,7)d at Tg = 0.1 
(upper panel) and Tg = 10 (lower panel). The solid line, 
for q = 1, corresponds to the usual Boltzmann distribu- 
tion. Also shown are results for the non-extensive dis- 
tributions for q — 0.5 (dotted line) and q — 2 (dashed 
line). One observes that, as for the charged particles 
case, with q < 1 higher kinetic energies are more proba- 
ble than in the extensive case (q = 1). With q > 1 high 
energies are less accessible than in the extensive case, 
and there is a cutoff beyond which no kinetic energy is 
reached. A noticeable difference form the case of charged 
particles is the absence of the Coulomb barrier and a cor- 
respondingly lack of suppression of the reaction rates at 
low energies. As the temperature increases, the relative 
difference between the Maxwell-Boltzmann and the non- 
Maxwellian distributions decrease appreciably. This will 
lead to a rather distinctive pattern of the reaction rates 
for charged compared to neutron induced reactions. 



For neutron-induced reactions, a good accuracy 
(within 0.1%) for the numerical calculation of the reac- 
tion rates with Eq. ( 12 ) is reached using the integration 
limits between E = and E — 20k bT. As an example 
we will now consider the reaction 7 Be(n,p) 7 Li. The en- 
ergy dependence of R{E) = S{E)\E for this reaction is 



non-extensive parameter q. 



shown in figure [8] The experimental data were collected 
from Refs. [5SH52]. 

Using the chi-square fit with a sum of polynomials and 
Breit-Wigners obtained to reproduce the data in figure [8j 
we show in figure[9]the reaction rates for 7 Bc(n,p) 7 Li as a 
function of the temperature Tg for different values of the 
non-extensive parameter q. The rates are given in terms 
of the natural logarithm of Na{o"v) (in units of cm 3 mol -1 
s^ 1 ). The solid curve corresponds to the usual Boltz- 
mann distribution, i.e., q = 1. The dashed and dotted 
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Figure 9: Reaction rates for 7 Be(n,p) 7 Li as a function of the 
temperature Tg for two different values of the non-extensive 
parameter q. The rates are given in terms of the logarithm 
of Na(cv) (in units of cm 3 moP 1 s _1 ). Results with the use 
of non-extensive distributions for q = 0.5 (dotted line) and 
q = 2 (dashed line) are shown. 
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Figure 10: Deuterium abundance. The solid curve is the re- 
sult obtained with the standard Maxwell distributions for the 
reaction rates. Results with the use of non-extensive distri- 
butions for q = 0.5 (dotted line) and q = 2 (dashed line) are 
shown. 



curves are obtained for q = 2 and q = 0.5, respectively. 
In contrast to reactions induced by charged particles, we 
now see strong deviations from the Maxwellian rate both 
for q > 1 and q < 1. For q < 1 the deviations are 
larger at small temperatures and decrease as the energy 
increase, tending asymptotically to the Maxwellian rate 
at large temperatures. This behavior can be understood 
from figure [n (for 7 Be(n,p) 7 Li the results are nearly the 
same as in Fig. ffh. At small temperatures, e.g. Tg =0.1, 
the distribution for q — 0.5 is strongly enhanced at large 
energies and the tendency is that the reaction rates in- 
crease at low temperatures. This enhancement disap- 
pears as the temperature increase (lower panel of figure 
[7]). For q = 2 the reaction rate is suppressed, although 
not as much as for the charged-induced reactions, the 
reason being due a compensation by an increase because 
of normalization at low energies. 

Having discussed the dependence of the reaction rates 
on the non-extensive parameter q for a few standard re- 
actions, we now consider the implications of the non- 
extensive statistics to the predictions of the BBN. It is 
clear from the results presented above that an apprecia- 
ble impact on the abundances of light elements will arise. 



calculated with a modified version of the standard BBN 
code derived from Refs. [63H65j. 

Although BBN nucleosynthesis can involve reac- 
tions up to the CNO cycle |B5], the most impor- 
tant reactions which can significantly affect the pre- 
dictions of the abundances of the light elements 
[ 4 Hc, D, 3 He, 7 Li] are n-decay, p(n,7)d, d(p,7) 3 He, 
d(d,n) 3 He, d(d,p)t, 3 He(n,p)t, t(d,n) 4 He, 3 He(d,p) 4 He, 
3 He(a,7) 7 Be, t(a, 7 ) 7 Li, 7 Be(n,p) 7 Li and 7 Li(p,a) 4 He. 
Except for these reactions, we have used the reaction 
rates needed for the remaining reactions from a compila- 
tion by NACRE [57] and that reported in Ref. [55]. For 
the 11 reactions mentioned above, we have collected data 
from Refs. 167-69], and references mentioned therein 
(data for n(p,7,d) reaction was taken from the on-line 
ENDF database [70] - see also HQ [7J), fitted the S- 
factors with a sum of polynomials and Breit-Wigner func- 
tions and calculated the reaction rates for Maxwellian 
and non-Maxwellian distributions. 



A. Elemental abundances 



IV. BBN WITH NON-EXTENSIVE STATISTICS 

The BBN is sensitive to certain parameters, including 
the baryon-to-photon ratio, number of neutrino families, 
and the neutron decay lifetime. We use the values rj — 
6.19 x 10~ 10 , N v = 3, and r„ = 878.5 s for the baryon- 
photon ratio, number of neutrino families, and neutron- 
day lifetime, respectively. Our BBN abundances were 



In figure [To] we show the calculated deuterium abun- 
dance. The solid curve is the result with the stan- 
dard Maxwell distributions for the reaction rates. Using 
the non-extensive distributions yields the dotted line for 
q = 0.5 and the dashed line for q = 2. It is interest- 
ing to observe that the deuterium abundances are only 
moderately modified due to the use of the non-extensive 
statistics for q = 0.5. Up to temperatures of the order 
of Tg ~ 1, the abundance for D/H tends to agree for the 
extensive and non-extensive statistics. This is due to the 
fact that any deuterium that is formed is immediately 
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destroyed (a situation known as the deuterium bottle- 
neck). But, as the temperature decreases, the reaction 
rates for the p(n,7)d reaction are considerably enhanced 
for q = 2 (see figure [7]) , and perhaps more importantly, 
they are strongly suppressed for all other reactions in- 
volving deuterium destruction, as clearly seen in figure [4] 
This creates an unexpected over abundance of deuterons 
for the non-extensive statistics with q = 2. The deu- 
terium, a very fragile isotope, is easily destroyed after 
the BBN and astrated. Its primordial abundance is de- 
termined from observations of interstellar clouds at high 
redshift, on the line of sight of distant quasars. These 
observations are scarce but allow to set an average value 
of D/H = 2.82±g;?g x 10" 5 [73]. The predictions for the 
D/H ratio with the q = 2 statistics (D/H = 5.70 x 10~ 3 ) 
is about a factor 200 larger than those from the standard 
BBN model, clearly in disagreement with the observa- 
tion. 
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Table I: Predictions of the BBN (with tjwmap = 6.2 x 10" 10 ) 
with Maxwellian and non-Maxwellian distributions compared 
with observations. All numbers have the same power of ten 
as in the last column. 





Maxwell Non-Max. 


Non-Max. 


Observation 




BBN 


q = 0.5 


g = 2 




4 He/H 


0.249 


0.243 


0.141 


0.2561 ±0.0108 


D/H 


2.62 


3.31 


570 


2.82l^(xHT 5 ) 


3 He/H 


0.98 


0.91 


69.1 


(l.l±0.2)(xl0 -5 ) 


7 Li/H 


4.39 


6.89 


356. 


(1.58±0.31)(xl0~ 10 ) 



A much more stringent constraint for elemental abun- 
dances is given by 4 He, which observations set at about 
4 He/H eF p = 0.2561 ± 0.0108 ^MM- The 4He abi- 
dance generated from our BBN calculation is plotted in 
The solid curve is the result obtained with 



figure 11 



the standard Maxwell distribution for the reaction rates. 
Using the non-extensive distributions yields the dotted 
line for q = 0.5 and the dashed line for q = 2. Again, 
the predicted abundances for q = 2 deviate substantially 
from standard BBN results. This time only about half of 
4 He is produced with the use of a non-extensive statistics 
with q = 2. The reason for this is the suppression of the 
reaction rates for formation of 4 He with q — 2 through 
the charged particle reactions t(d,n) 4 He, 3 He(d,p) 4 He. 



Figure 11: 4 He abundance. The solid curve is the result ob- 
tained with the standard Maxwell distributions for the reac- 
tion rates. Results with the use of non-extensive distributions 
for q — 0.5 (dotted line) and q — 2 (dashed line) are also 
shown. 
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Figure 12: 3 He abundance. The solid curve is the result ob- 
tained with the standard Maxwell distributions for the reac- 
tion rates. Results with the use of non-extensive distributions 
for q — 0.5 (dotted line) and q — 2 (dashed line) are also 
shown. 



A strong impact of using non-extensive statistics for 
both q = 0.5 and q — 2 values of the non-extensive 

for the 3 He abundance. 
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parameter is seen in figure 
While for q = 2 there is an overshooting in the produc- 
tion of 3 He, for q = 0.5 one finds a smaller value than the 
one predicted by the standard BBN. This is due to the 
distinct results for the destruction of 3 He through the 
reaction 3 He(n,p)t, which is enhanced for q = 0.5 and 
suppressed for q = 2, in the same way as it happens for 
the reaction 7 Be(n,p) 7 Li, shown in figure [8] 3 He is both 



produced and destroyed in stars and its abundance is still 
subject to large uncertainties, 3 He/H = (1.1 ±0.2) x 10~ 5 

E7l ESI. 



Non-extensive statistics for both q = 0.5 and q — 2 val- 
ues also alter substantially the 7 Li abundance, as shown 



in figure 13 For both values of the non-extensive param- 
eter q — 2, and q = 0.5, there is an overshooting in the 
production of 7 Li. The increase in the 7 Li abundance 
is more accentuated for q = 2. The lithium problem 
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is associated with a smaller value of the observed 7 Li 
abundance as compared the predictions of BBN. There 
has been many attempts to solve this problem by test- 
ing all kinds of modifications of the parameters of the 
BBN or the physics behind it (a sample of this liter- 
ature is found in Refs. E2J EU 111] and refer- 
ences therein). In the present case, the use of a non- 
Maxwellian velocity distribution seems to worsen this 
scenario. A recent analysis yields the observational value 
of 7 Li/H = (1.58 ± 0.31) x 1(T 10 [85]. 



B. Sensitivity study 

We have calculated a window of opportunity for the 
non-extensive parameter q with which one can reproduce 
the observed abundance of light elements. We chose the 
data for the abundances, Yi, of 4 He/H, D/H, 3 He/H and 
7 Li as reference. We then applied the ordinary x 2 statis- 
tics, defined by the minimization of 



4^ CT,; 



(28) 



where Yi(q) are the abundances obtained with the non- 
extensive statistics with parameter q, Yi(obs) are the ob- 
served abundances, and (Ji the errors for each datum, and 
the sum is over all data mentioned in Table HI From this 
chi-square fit we conclude that q = l^o'i2 ^ s compatible 
with observations. 

No attempt has been made to determine which ele- 
ment dominates the constraint on q. This might be im- 
portant for a detailed study of the elemental abundance 
influence from nuclear physics inputs, namely, the uncer- 
tainty of the reaction cross sections. A study along these 
lines might be carried in a similar fashion as described in 
Ref. [75] . Weights on the reliability of observational data 
should also be considered for a more detailed analysis. 
For example, constraints arising from 3 He/H abundance 
may not be considered trustworthy because of uncertain 
galactic chemical evolution. On the other hand, a con- 
straint from the observation of 3 He/D is more robust, and 
so on. Based on our discussion in Section II, it is more im- 
portant to determine how the non-extensive statistics can 
modify more stringent conditions during the big bang, 
such as the modification of weak-decay rates and its in- 
fluence on the n/p ratio which strongly affects the 4 He 
abundance. 



V. CONCLUSIONS 

In table [I] we present results for the predictions of the 
BBN with Maxwellian and non-Maxwellian distributions. 
The predictions are compared with data from observa- 
tions reported in the literature. It is evident that the re- 
sults obtained with the non-extensive statistics strongly 
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Figure 13: 7 Li abundance. The solid curve is the result ob- 
tained with the standard Maxwell distributions for the reac- 
tion rates. Results with the use of non-extensive distributions 
for q — 0.5 (dotted line) and q — 2 (dashed line) are also 
shown. 



disagree with the data. The overabundance of 7 Li com- 
pared to observation gets worse if q > 1. The three 
light elements D, 4 He and 7 Li constrain the primordial 
abundances rather well. For all these abundances, a non- 
extensive statistics with q > 1 leads to a greater discrep- 
ancy with the experimental data. 

Except for the case of 3 He the use of non-extensive 
statistics with q < 0.5 does not rule out its validity 
when the non-Maxwellian BBN results are compared to 
observations. 3 He is at present only accessible in our 
Galaxy's interstellar medium. This means that it can- 
not be measured at low metallicity, a requirement to 
make a fair comparison to the primordial generation of 
light elements. This also means that the primordial 3 He 
abundance cannot be determined reliably. The result 
presented for the 3 He abundance in table [i] is quoted 
from Ref. [77]. Notice that our analysis does not in- 
clude the changes that the non-extensive statistics would 
bring to the n/p conversion rates. The electron distri- 
butions would also be expected to change accordingly. 
This would change the freeze out temperature and a cor- 
responding influence on the 4 He abundance. 



We conclude that it does not seem possible to change 
the Maxwell-Boltzmann statistics to reproduce the ob- 
served abundance of light elements in the universe with- 
out destroying many other successful predictions of big 
bang nucleosynthesis. A chi-square fit of our calcula- 
tions with the observations of elemental abundance con- 
cludes that the non-extensive parameter is constrained to 
q = lio 12' This means that, should a non-Maxwellian 
distribution due to the use of the Tsallis non-extensive 
statistics be confirmed (with a sizable deviation from 
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q = 1), our understanding of the cosmic evolution of the 
universe would have to be significantly changed. 
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